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Abstract 

Dissipation of magnetic fields in Hall plasma of neutron star crusts 
may power persistent high energy emission of a class of strongly mag- 
netized neutrons stars, magnetars. We consider development of a dis- 
sipative tearing mode in Ifall plasma (electron MHD) and find that its 
growth rate increases with the wave number of perturbations, reaching 
a maximum value intermediate between resistive and Hall times th, 
r ~ l/s/TrTH- We argue that the tearing mode may be the principal 
mechanism by which strong magnetic fields are dissipated in magnetars 
on times scale of ~ 10'* — 10^ yrs powering the persistent X-ray emisison. 

1 Introduction 

Evolution of neutron stars ' magnetic fields is one of the principal issues in 
neutron star research. On the one hand, observationally, there is no evidence 
that fields decay on long time scales, of the order of million years [1]. Yet, this 
constrains the dipolar component of a magnetic field, generated, presumably, by 
currents flowing in the superconductive core. The fate of a magnetic field in the 
crust may be different. It is tempting to relate an X-ray activity of magnetars, 
strongly magnetized neutron stars , to the decay of a crustal magnetic field 
created by a dynamo action during the neutron star birth [2, 3]. 

One of the central problems is how fast a magnetic field can be dissipated. 
Conventionally, the magnetic field decay time scale is assumed to be of the order 
of the resistive time scale, given by the size of the system L (of the current- 
carrying layer, to be more precise) and resistivity 77, L'^ /v- In fact, plasma 
has a number of ways to dissipate magnetic fleld on time scales much shorter 
than Tr (otherwise we would never see a Solar flare, since the resistive time scale 
in the Solar corona is longer than the age of the Universe). Typically, fast mag- 
netic fleld dissipation is achieved through formation of small scale structures, 
with correspondingly short dissipation time scales. One possibility to create 
small scale structures is through non-linear interaction of waves, either locally 
(in phase space) through formation of a Hall cascade [4] , or non- locally due to 
large scale motions leading to "wave overturn" [5, 6, 7, 8, 9]. 
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Here we discuss another way to form small scale dissipative current layers 
via development of a tearing mode in Hall plasma. In Hall plasma (sometimes 
called electron MHD, EMHD below [5]) ions are assumed to be motionless, 
providing a neutralizing background for electron fluid. Tearing mode [10] is 
the principal resistive instability of current-currying plasma and is one of the 
key factors leading to an explosive release of magnetic energy in Solar flares 
[11] and magneto-tail [12]. It is also important for TOKAMAK discharges like 
sawtooth oscillations and major disruptions [13]. Tearing mode in regime is 
an important ingredient of modern reconnection models [7, 14]. During the 
development of the tearing mode, perturbations of a plasma equilibrium lead to 
formation of current sheets where both the time scales for diffusion may be short 
and, in addition, resistivity may be enhanced due to the development of plasma 
turbulence (anomalous resistivity) . A current layer tends to be locally unstable 
to transverse, k • Bq = 0, perturbations (k is the wave vector of perturbations 
and Bo is initial magnetic field). Qualitatively, divergence at points k • Bq = 
is related to the fact that Alfven velocity becomes zero at these points, so that 
perturbations are effectively piling up. Since whistler waves, normal modes in 
Hall plasma, also have phase speed equal to zero at points where k • Bq = 0, we 
might expect a somewhat similar behavior. This is indeed what we show below. 

2 Tearing mode in Hall plasma 
2.1 Teeiring due to electron resistivity 

In deriving growth rates for the tearing mode in Hall plasma we follow a standard 
scheme of electron-ion plasma [15]. On scales much larger that electron skin 
depth, the equation of resistive EMHD is (for example [5]) 



where rj is resistivity, n is plasma density and other notations are standard. In 
what follows we assume density to be constant. 

We assume that initial magnetic field configuration is stable both on Alfven 
time scale (implying a dynamical stability) and on Hall time scale. In doing so 
we wish to single out the effects of the tearing mode as opposed to dynamical 
instabilities. Hall cascade etc. 

Let the initial magnetic field correspond to current layer B = Bo{fx{z)ex + 
fy{z)ey) with /a; ^ ±1 at 2 ^ ±00 and /j; = in the center of the layer z = 0. 
(If the layer is force-free, then f^+fy = 1; this is not required for tearing 
instability). Next we impose perturbation of the field in the form 



where ^ cx e ''"^ '^"^^ is vector potential. We foresee instability at points 
k • Bq = 0, so to simplify mathematics, we assume ky = 0, so that a current 




(1) 



(5B = V X ^ 



(2) 
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sheet will be created at ^; = 0. The equations for evolution of perturbations are 



< + i'^MU^ + M'J - ^yf'^ - CU - ILC - fyO + 

^{ik.^. - O + {klUy + - Uy) 

-V {kM + O + + ^{klC - e^)) = 

-ik^cj^y + {k.M. + iUy + iUC) + Vk. {k% - O = (3) 

The two scales of interest that appear in the equations are 




= {Lk^dy-^ (4) 

C^kxCOB tOw 

where ojyj = {ckx)'^UB/i^p is the frequency of whistler waves. Generally, dis- 
persion relation of whistler waves ui^ = (^kkxOJs/'jJp = c^fc(k • wb)/^^, where 
k = i/fc|"+fcf . In an inhomogeneous plasma of our choice Wu, = at a point 
2 = 0. For comparison, in c-i plasma 5* = \J L&\lij\I uj a, with uja = kxVA- In the 
above equations |a;| is an absolution value of the (generally complex) frequency 
uj. As it turns out, real part of to is not important (see below), so |w| ^ F, where 
r is a growth rate of instability. 

A standard method in describing the evolution of tearing mode is similar 
to the boundary layer problem. It involves separation of a current layer into 
"bulk", where derivatives arc small and resistivity is not important, and a nar- 
row "boundary layer" where derivatives and resistivity may be large. Two dif- 
ferent approximations are done in each layer - ideal and weakly varying plasma 
in the bulk and a narrow resistive sublayer, with two solutions matched contin- 
uously. 

In the outside region of smooth field evolution, where we can neglect resis- 
tivity, the system (3) reduces to 

We are interested in the slow evolution of plasma and thus neglect terms propor- 
tional to iv. The structure of magnetic field in the ideal region is then described 
by 

= (kl + iy (6) 
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This clearly shows that instability is driven by term, the second derivative 
of magnetic field. 

The simplest choice of the background field in (6) is fx = smz/ L, z < L, in 
which case 



^y = {kl-^yy (7) 



with solution 



iy = Ci cos I ^ 1 + Ca sm I ^ 1 (8) 

Outside of the current sheet, z > L, = e~'^»'^ and the matching gives 

kxL sin a/1 — fc^L^ 



Ci = e-''"^ cos - klL-^ 



Thus, magnetic field = kx^y is continuous at ^; = 0, but its derivative 
experiences a jump 

A=[lnB,]--^ (10) 

For instability it is required that A > 0. (Recall, for comparison, that in 
electron-ion plasma A l/{kL?)). 

In the resistive sub-layer, z ^ L, we have /a, — > , — > 0, so that the z 
component of (3) gives 

iy - -^^y (11) 

where S"^ = tj/T, T = Im(w). In the limit kxS <C 1 this gives 

^y= cosh z/S (12) 

where a boundary condition = at 2; = was used. (Structure of the resistive 
sub-layer is the same as in electron-ion plasma). The inner and outer solutions 
match at 

S* = AS"^ (13) 

where 6* is the thickness of the resistive sub-layer. 
Prom Eqns (4), (13) we get the growth rate 



^ _ c^kxS'^A'^UB _ f (^kxJpJB\ 



(14) 



This is consistent with Eq. 50 in Ref. [7], and Eq. (68) in Ref. [16]; see also 

[17, 18]. If we define phase velocity of whistler modes as = kxC^cos/i^p, then 
Eq. (14) takes the form 

r = ^ (15) 
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where resistive time scale \STr = L"^ /rj and whistler time scale = L/vw 

1 /2 

Growth rate (14) increases as K'\ This is in stark contrast to electron-ion 
plasma, where growth rate decreases with kx (as oc k~'^/^, e.g. [15] Eq. 7.123). 
Thus, in Hall plasma the maximum growth rate is reached at kx ~ 1/i^, and 
becomes 

r = ^ = ^ (16) 

where we defined Hall time 

7" 2, ,2 

= ^ (17) 

Relation (14) gives the growth rate of tearing mode in Hall plasma. Thus, 
the time scale for development of a tearing mode is intermediate between the 
resistive and Hall time scales, and thus is much shorter than resistive time 
scale. This expression reminds of the maximum growth rate of tearing mode 
in electron-ion plasma, which is intermediate between resistive and Alfven time 
scales. 

Growth rate (16) may also be written as 

( ^ V '^-B 



\ r } r- (18) 

where we introduced an effective Lundquist number associated with Hall time 

S = vhL/t], vh = L/th 

2.2 Teeiring due to inertial resistivity 

In a highly conducting plasma, electron inertia may play a role of resistivity, 
providing a relation between the electric field and the current. It is this regime 
that most works on tearing mode in Hall plasma addressed so far [7]. Since in 
this case the width of the current sheet may become comparable to electron skin 
depth, one has to take into account additional terms in the Ohm's law and it's 
consequence, Eq. (1) If we neglect this possibility (this is justified since we get 
correct expression for the growth rate, sec below), the effective resistivity in the 
coUisionless regime becomes rjefj = r(c/ci;p) [19, 7]. Thus, inertial resistivity 
dominates when T{c/(jjpY > t], which can be written as {cj Llo^Y Tr j t^. 
The corresponding tearing mode growth rate is 

\ 4 4, 



This is the same as the growth for electron tearing mode cited in Ref. [20] . The 
maximum rate is reached at fc ~ l/L 

^'""^ = (z^) 

The ratio of resistive (16) and inertial (20) maximum growth rates is {c/ LujpY^/S. 
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3 Application to neutron stars 



High energy emission of strongly magnetized neutron stars (magnetars) is pow- 
ered by dissipation of magnetic field [3, 21, 22]. Strong magnetic fields, of the 
order of 10^^ G, may be created by a dynamo mechanism, e.g. of the a — Q, type, 
operating at birth of neutron stars [2]. After ~ 100 sees from the birth, the crust 
of the neutron star solidifies. This time is much longer than Alfvcn crossing time 
'--^ 0.1 — 1 s, so that after the end of the turbulent motion and before the crust 
formation, magnetic field in a star should evolve to some minimum energy state 
allowed by the system. 

After the crust solidifies, the evolution of magnetic field will proceed due to 
non-dissipative Hall effect, Ohmic resistivity, and, as we argue in this paper, due 
to development of tearing instability. In order to estimate the corresponding 
times scales, we use resistivity r] = /A-ku with conductivity a ~ lO^^p^j'^ [23] 
evaluated at the neutron drip point p 10^^ g/cm~^. The Hall th, resistive 
and tearing time scales (given by inverse of the tearing mode growth rate (14)) 
become 

TH = Ax 10^LIB^Ipi2 yrs 
Tr = 4:X IG'^i^pJf yrs 

l/T = 4xlO'Llp%'B-^'/^yvs (21) 

where a standard subscript notation, e.g. L4 = (L/lC^cm), has been adopted. 
This indicates that growth rate of the tearing mode in Hall plasma is of the 
same order as the activity time of magnetars. Also, for higher magnetic field 
the instability time is shorter while corresponding magnetic energy is larger, 
consistent with the fact that only high field neutron stars exhibit magnetar-like 
activity. One can also check that in neutron star crusts resistive effects well 
dominate over inertial ones. 

We envision a magnetar scenario similar to the ones proposed in Refs [24, 25] . 
After the turbulence seizes, initially stable magnetic field configuration forms 
which is then dissipated due to development of tearing mode on time scale 
~ 10** — 10^ yrs. This leads to Lorentz force disbalance in the crust, which 
initially can be compensated by crust tensile strength, but eventually is released 
through crust deformation leading to bursts and flares. (Release of crustal stress 
can be either through cracks [2] or through plastic deformations [26]). 

Since small scale currents are dissipated on shorter time scales, while larger 
currents take longer to dissipate, we expect that giant flares are more common 
in older magnetars. (This is also a natural consequence of the torus formation 
[27], in which case small scale magnetic fields are dissipated soon after the birth 
of a neutron star ) . This may provide a resolution to the energy budget problem 
for magnetar giant flares: the giant flare of SGR 1806 - 20 emitted 2 x lO"'^ 
ergs in high energy 7-rays [28] , and associated mechanical energy is expected to 
be even larger. If a typical giant fiare recurrence time is ^ 30 yrs (the period 
of active monitoring of high energy sky) and activity time is several thousand 
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years (characteristic age) the magnetic field, which powers the flares, has to be 
larger than 3 x 10^^ G, which is on the verge of being uncomfortably large. 

4 Conclusion 

In this paper we first considered development of a resistive tearing mode in Hall 
plasma. Tearing mode develops on time scale intermediate between Hall and 
resistive time scales. Qualitatively, this might have been expected in analogy 
with the tearing mode in electro-ion plasma. Unlike the case of electron-ion 
plasma, the growth rate (14) increases with k, reaching maximum for k 1/L, 
the thickness of the current layer. This is due, qualitatively, to the dispersion 
of whistler modes, for which phase velocity increases with k. 

We then argue that the development of the tearing mode in the crust of 
strongly magnetized neutron stars, magnetars, powers their high energy emis- 
sion on time scales ~ 10'' — 10^ years. One of the principal uncertainty is the 
length scale L of the current sheets created by the dynamo, as all the time scales 
are strongly dependent on it, Eq. (21). For our choice of parameters the Hall 
time scale is th ~ lO"' yrs. According to Ref. [23], this, in fact, applies to 
a very broad density range, p ~ 10^^ — 10^^ g cm~^. Ohmic decay times are 
strongly dependent on the temperature, with the value that we used, tr ~ 10^ 
yrs., giving a reasonable approximation. 

Overall, the evolution of a crustal magnetic field is bound to be a combi- 
nation of a dissipativc tearing mode and a nondissipativc Hall cascade. The 
Hall turbulent cascade can create small scale structures which in turn become 
dissipative [4] . The time it takes for a cascade to propagate down to dissipative 
scale is, typically the large eddy overturn time times the logarithm of the ratio 
of outer to inner scales, ~ Ti/ In ^ th [29] , where is the dissipa- 

tive (resistive) scale. As the ratio L/Lin is large, it takes many Hall times to 
dissipate magnetic energy. Numerical simulations [30] indeed seem to indicate 
that the transfer of energy to the higher harmonics is not sufficient to accelerate 
significantly the decay of the original field. 

Relative importance of Hall cascade and tearing mode depend on details of 
the structure of magnetic field before solidification of the crust. This is yet an 
unresolved issue, as is exemplified by a long standing problem of stability of 
magnetic field in stars [31, 32, 27, 33, 34]. It appears that evolution at interme- 
diate times (longer than Alfven crossing time but shorter than dissipative times) 
depends on magnetic helicity, either at the end of the turbulent phase or due 
to helicity imbalance due to loss through the surface). For certain parameters, 
internal magnetic field relaxes to a complicated torus-like form with comparable 
toroidal and poloidal magnetic fields. 

Our results provides an elegant explanation to observation in Ref. [35] (see 
also [36]) where it is found that formation of dissipative current sheets requires 
at least quadratic dependence of the background field on the position. Eq. 
(6) clearly shows that instability is driven by the second spacial derivative of 
the field. Thus, what is calculated in Refs [35] is a tearing mode. Finally, in 
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Rcf. [8] a somewhat related problem was considered, the formation of a current 
sheet in presence of a steep density gradient. In contrast, the tearing mode 
considered in this paper does not require density gradient, all that is needed 
is an inhomogencous magnetic field. In addition, limitation of purely toroidal 
field, assumed in [8] seems to overestimate the dissipation [37]. 

Finally, the model presented here may be directly probed in the laboratory, 
in particular at the field-reversed configuration experiment at UCLA [38] . Tear- 
ing and formation of magnetic islands, a typical consequence of the nonlinear 
development of the tearing mode, are clearly seen (e.g. Fig 9 in Ref [38]). 
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Holler bach, Ulrich Geppert, Matthias Reinhardt, Gennady Shvets and Dmitri 
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